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7 Abstract. Let P be a set of n pairwise disjoint unit disks in the plane. 

^ 8 We describe how to build a data structure for V so that for any point 

^ 9 set P containing exactly one point from each disk, we can quickly find 

10 the onion decomposition (convex layers) of P. 

, 11 Our data structure can be built in 0(n\og7i) time and has linear size. 

12 Given P, we can find its onion decomposition in 0{nlogk} time, where 

13 k is the number of layers. We also provide a matching lower bound. 

14 Our solution is based on a recursive space decomposition, combined with 

15 a fast algorithm to compute the union of two disjoint onion decomposi- 

16 tions. 



u 

17 1 Introduction 



18 Let P be a planar n-point set. Take the convex hull of P and remove it; repeat 

19 until P becomes empty. This process is called onion peeling, and the resulting 

20 decomposition of P into convex polygons is the onion decomposition, or onion for 
Cf^ 21 short, of P. It can be computed in 0(n log n) time [6]. Onions provide a natural, 
•/^ 22 more robust, generalization of the convex hull, and they have applications in 

23 pattern recognition, statistics, and planar halfspace range searching [7,14 22 



24 Recently, a new paradigm has emerged for modeling data imprecision. Sup- 

25 pose we need to compute some interesting property of a planar point set. Suppose 
I 26 further that we have some advance knowledge about the possible locations of the 

^> 27 points, e.g., from an imprecise sensor measurement. We would like to preprocess 

28 this information, so that once the precise inputs are available, we can obtain 

29 our structure faster. We will study the complexity of computing onions in this 

30 framework. 

31 1.1 Related Work 

32 The notion of onion layer decompositions first appears in the computational 

33 statistics literature [l4] , and several rather brute-force algorithms to compute it 

34 have been suggested (see [9] and the references therein). In the computational 

35 geometry community, Overmars and van Leeuwen [2l] presented the first near- 

36 linear time algorithm, requiring 0{nlog'^ n) time. Chazelle [H] improved this 



37 to an optimal 0(n log n) time algorithm. Nielsen 20 gave an output-sensitive 



(a) (b) 
Fig. 1. (a) Two disjoint onions, (b) Their union. 



1 algorithm to compute only the outermost k layers in 0{nloghk) time, where hk 

2 is the number of vertices participating on the outermost k layers. In M.^, Chan [sj 

3 described an 0{nlog^ n) expected time algorithm. 

4 The framework for preprocessing regions that represent points was first in- 



5 troduced by Held and Mitchell 12 , who show how to store a set of disjoint 

6 unit disks in a data structure such that any point set containing one point from 

7 each disk can be triangulated in linear time. This result was later extended to 



8 arbitrary disjoint regions in the plane by van Kreveld et al. 16 . LofHer and 

9 Snoeyink first showed that the Delaunay triangulation (or its dual, the Voronoi 

10 diagram) can also be computed in linear time after preprocessing a set of disjoint 

11 unit disks 17 . This result was later extended by Buchin et al. [2], and Devillers 

12 gives a practical alternative [sj. Ezra and Mulzer 10 show how to preprocess 

13 a set of lines in the plane such that the convex hull of a set of points with one 

14 point on each line can be computed faster than n log n time. 

15 These results also relate to the update complexity model. In this paradigm, 

16 the input values or points come with some uncertainty, but it is assumed that 

17 during the execution of the algorithm, the values or locations can be obtained 

18 exactly, or with increased precision, at a certain cost. The goal is then to compute 

19 a certain combinatorial property or structure of the precise set of points, while 

20 minimising the cost of the updates made by the algorithm [3| [Ilj[T3|[23] . 



21 1.2 Results 

22 We begin by showing that the union of two disjoint onions can be computed in 

23 0{n + k^ log ri) time, where k is the number of layers in the resulting onion. 

24 We apply this algorithm to obtain an efficient solution to the onion prepro- 

25 cessing problem mentioned in the introduction. Given n pairwise disjoint unit 

26 disks that model an imprecise point set, we build a data structure of size 0{n) 

27 such that the onion decomposition of an instance can be retrieved in 0{n\ogk) 

28 time, where k is the number of layers in the resulting onion. We present several 

29 preprocessing algorithms. The first is very simple and achieves O(nlogn) ex- 

30 pected time. The second and third algorithm make this guarantee deterministic, 

31 at the cost of worse constants and/or a more involved algorithm. 

32 We also show that the dependence on k is necessary: in the worst case, any 

33 algorithm can be forced to take Q{n\ogk) time on some instances. 
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1 2 Preliminaries and Definitions 

2 Let P be a set of n points in M^. The onion decomposition, or onion, of P, is 

3 the sequence ^(P) of pairwise disjoint convex polygons with vertices from P, 

4 constructed recursively as follows: if P 7^ 0, we set #(P) := {ch(P)} U #i(P \ 

5 ch(P)), where ch(P) is the convex hull of P; if P = 0, then #(P) := An 

6 element of i#(P) is called a layer of P. 

7 Let 2? be a set of disjoint unit disks in M.~. We say a point set P is a sample 

8 from V if every disk in V contains exactly one point from P. We write log for 

9 the logarithm with base 2. 



10 3 The Algorithm 



11 Our algorithm requires several pieces, to be described in the following sections. 



12 3.1 Unions of Onions 



13 Suppose we have two point sets P and Q, together with their onions. We show 

14 how to find #>(P U Q) quickly, given that ^(P) and ^{Q) are disjoint. Deleting 

15 points can only decrease the number of layers, so: 

16 Observation 3.1 Let P,Q Q M? . Then ^{P) and ^{Q) cannot have more lay- 

17 ers than ^{P U Q) . □ 

18 The following lemma constitutes the main ingredient of our onion-union al- 

19 gorithm. A convex chain is any connected subset of a convex closed curve. 

20 Lemma 3.1. Let A and B be two non-crossing convex chains. We can find 

21 ch{AlJ B) in 0{\og n) time. 



22 Proof. Since A and B do not cross, the pieces of A and B that appear on 

23 ch(A U B) are both connected: otherwise, ch(yl U B) would contain four points 

24 belonging to A, B, A, and B, in that order. However, the points on A must be 

25 connected inside c\i{A U B); as do the points on B. Thus, the chains A and B 

26 cross, which is impossible. Since A and B are convex chains, we can compute 

27 ch(A),ch(P) in O(logTi) time. Furthermore, since A and B are disjoint, we can 

28 also, in O(logn) time, make sure that ch(A) n ch(P) = 0, by removing parts 

29 from A or B, if necessary. Now we can find the bitangents of ch(A) and ch.{B) 



30 in logarithmic time 15 . □ 



31 Lemma 3.2. Suppose #i(P) has k layers. Let A be the outer layer of <$( 

32 andp.q be two vertices of A. Let Ai be the points on A between p and q, going 

33 counter-clockwise. We can find %{P\ Ai) in 0{k log n) time. 
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(a) 



(b) 



Fig. 2. (a) A half-eaten onion; (b) the restored onion. 
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Proof. The points p and q partition A into two pieces, Ai and A2. Let B be the 
second layer of ^(P). The outer layer of ^(P\Ai) is the convex hull of P\Ai, 



i.e., the convex hull of .42 and B. By Lemma 3.1 we can find it in 0(log n) time. 
Let p'.q' E P he the points on B where the outer layer of ^(P \ Ai) connects. 
We remove the part between p' and q' from B, and use recursion to compute 
the remaining layers of ^{P \ Ai) in 0{{k — 1) logn) time; see Figure[2] □ 



7 We conclude with the main theorem of this section: 



Theorem 3.2. Let P and Q be two planar point sets of total size n. Suppose 
that and #>((5) are disjoint. We can find the onion #(PU Q) in 0{k'^ log n) 

time, where k is the resulting number of layers. 



Proof. By Observation 



3.1 



(P) and 



each have at most k layers. We use 



Lemma 3.1 to find cli(PU Q) in 0(log n) time. By Lemma 3.2 the remainders of 
^(P) and <^{Q) can be restored to proper onions in 0(A;logn) time. The result 
follows by induction. □ 



15 3.2 Space Decomposition Trees 

16 We now describe how to preprocess the disks in V for fast divide-and-conquer. 

17 A space decomposition tree (SDT) T is a rooted binary tree where each node 

18 V is associated with a planar region R.^,. The root corresponds to all of M'^; for 

19 each leaf v of T, the region R^, intersects only a constant number of disks in 

20 T). Furthermore, each inner node v in T is associated with a directed line so 

21 that if u is the left child and w the right child of w, then i?„ := J?„ n and 

22 Ru, := Rv n £^ . Here, i'^ is the halfplane to the left of C^. and the halfplane 

23 to the right of £.„ ; see Figure [3] 

24 Let a, /3 e (0, 1), and let T be an SDT. For a node v of T, let f/„ denote the 

25 number of disks in T> that intersect Ry. We call T an (q;,/3)-SDT for T> if for 

26 every inner node v we have that (i) the line intersects at most d^l disks in 

27 and (ii) d„, d^ < ady, where u and w are the children of v. 
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1 Lemma 3.3. Let T be an [a, l3)-SDT. The tree T has height O(logn) and 0{n) 

2 nodes. Furthermore, X^ugT^" ~ 0{nlogn). 

3 Proof. The fact that T has height O(logn) is immediate from property (ii) of 

4 an (a,/3)-SDT. For i = 0, . . . ,logn, let V, := {v e T \ d., e [2',2'+i)}, the set 

5 of nodes whose regions intersect between 2' and 2*+^ disks. Note that the sets 

6 Vi constitute a partition of the nodes. Let Vi C Vi be the nodes in Vi whose 

7 parent is not in Vi. By property (ii) again, the c?„ along any root-leaf path in T 

8 are monotonically decreasing, so the nodes in Vi are unrelated (i.e., no node in 

9 Vi is an ancestor or descendant of another node in V). Furthermore, the nodes 

10 in Vi induce in T a forest F, such that each tree in Fi has a root from Vi and 

11 constant height (depending on q). 

12 Let Di :— J2vev- claim that for i = 0, . . . , log n, we have 

log n 

D,<nY[{l + c2^^^-'^), (1) 

13 for some large enough constant c. Indeed, consider a node v G Vj. As noted 

14 above, v is the root of a tree Fy of constant height in the forest induced by Vj. 

15 By property (i), any node u in this subtree adds at most < 2^^+^''^ additional 

16 disk intersections (i.e., da + di, < du + 2^^^^''^ , where a, b are the children of 

17 u). Since Fy has constant size, the total increase in disk intersections in F„ is 

18 thus at most c'2'^+^'^, for some constant c'. Since dy > 2^, it follows that the 

19 number of disk intersections increases multiplicatively by a factor of at most 

20 1 + c'2(J+i)'3/2-' < 1 + c2J('^-i', for some constant c. The trees Fy partition T 

21 and the root intersects n disks, so for the nodes in V., the total number of disk 

22 intersections has increased by a factor of at most nj-^i"(l + c2-'^^~^^), giving 

23 ([T]). The product in ([T]) is easily estimated: 

log n 

A < n n (1 + c2^(^-i^) < ne^i- = ne^^^) = 0(n), 

24 since /3 < 1. Hence, each set Vi has at most 0(n/2') nodes for i = 1, . . . , logn. 

25 The total size of all Vi is 0(n). Since each v E Vi lies in a constant size subtree 

26 rooted at a w G Vi, it follows that T has 0{n) nodes. For the same reason, we 

27 get that X^dst'^i' ~ 0{nlogn}. □ 

28 Now there are several ways to obtain an (a,/3)-SDT for V. A very simple 

29 construction is based on the following lemma, which is an algorithmic version of 

30 a result by Alon et al. |2| Theorem 1.2]. See Section |4] for alternative approaches. 

31 Lemma 3.4. There exists a constant c > 0, so that for any setT> of m congruent 

32 nonoverlapping disks in the plane, there is a line £ with at least m/2 — c^/m log m 

33 disks completely to each side of it. We can find £ in 0{ni) expected time. 
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Fig. 3. A space decomposition tree for 21 unit disks. 



1 Proof. Our proof closely follows Alon et al. [2j Section 2]. Set := m/ log m\ , 

2 and pick a random integer z between 1 and r/2. Find a line (. whose angle with 

3 the .r-axis is [z/r]?: and that has \m/2\ disk centers on each side. Given r, 

4 we can find C in 0{m) time by a median computation. The proof by Alon et al. 

5 implies that with probability at least 1/2 over the choice of the line (' intersects 

6 at most c^m log m disks in V, for some constant c > 0. Thus, we need two tries 

7 in expectation to find a good line f. The expected running time is 0{m). □ 



8 To obtain a (1/2 + e, 1/2 + £)-SDT T for V, we apply Lemma |3 . 4| recursively 

9 until the region for each node intersects only a constant number of disks. Since 
10 the expected running time per node is linear in the number of intersected disks, 



11 Lemma 3.3 shows that the total expected running time is 0(ri log n). 



12 By Lemma 3.3 the leaves of T induce a planar subdivision Gj- with 0{n) 

13 faces. We add a large enough bounding box to Gt and triangulate the resulting 

14 graph. Since Gt is planar, the triangulation has complexity 0{n) and can be 

15 computed in the same time (no need for heavy machinery — all faces of Gt are 

16 convex). With each disk in V, we store the list of triangles that intersect it (recall 

17 that each triangle intersects at most one disk). This again takes 0{n) time and 

18 space. We conclude with the main theorem of this section: 

Theorem 3.3. Let V be a set of a disjoint unit disks in Ml^ . In O(nlogn) 
expected time, we can construct an (1/2 + e, 1/2 + e) space partition tree T for 
V. Furthermore, for each disk D £ T) , we have a list of triangles Td that cover 
the leaf regions of T that intersect D. □ 



19 3.3 Processing a Precise Input 

20 Suppose we have an (q, /:f)-SDT together with a point location structure as 

21 in Theorem 3.3 Let P be a sample from V. Suppose first that we know k, 



22 the number of layers in m{P). For each input point pi, let Di <E V he the 

23 corresponding disk. We check all triangles in Tq. , until we find the one that 
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1 contains pi. Since there are 0{n) triangles, this takes 0{n) time. Afterwards, we 

2 know for each point in P the leaf of T that contains it. 

3 For each node v of T, let n,^ be the number of points in the subtree rooted 

4 at V. We can compute the n.^^s in total time 0(ti) by a postorder traversal of 

5 T. The upper tree of T consists of all nodes v with n„ > k^. Each leaf of 

6 Tu corresponds to a subset of P with 0{k~) points. For each such subset, we 

7 use Chazelle's algorithm [6| to find its onion decomposition in O(fc^logfc) time. 

8 Since the subsets are disjoint, this takes O(nlogfc) total time. Now, in o rder to 

9 obtain ^(P), we perform a postorder traversal of T^, using Theorem 3.2 in each 

10 node to unite the onions of its children. This gives ^{P} at the root. 

11 The time for the onion union at a node v is 0{k^ logn^,). We claim that for 

12 i ~ 2\ogk, . . . , log n, the upper tree T,, contains at most 0{n/T-) nodes v with 

13 riy G [2*, 2*+^). Given the claim, the total work is proportional to 

log n log n . 

Efc'log""< E J^k\^ + l)^nk^ £ l±1^0(„logfc), 

■ueTu 1=2 log fc i=21ogfc 

14 since the series Y^i^2\agki'^ + is dominated by the first term (logfc)/fc^. 

15 It remains to prove the claim. Fix i £ {2 log fc, . . . , log n} and let Vi be the 

16 nodes in Ta with E [2;, 2*+^), whose parents have Uy > 2*+^. Since the 

17 nodes in Vi represent disjoint subsets of P, we have \Vi,\ < n/2*. Furthermore, 

18 by property (i) of an (q:,/3)-SDT , both children u'i, W2 for every node v G Tu 

19 have 71^,^,11^., < onii,, so that after 0(1) levels, all descendants w of n G V have 

20 < 2'. The claim follows. 

21 So far, we have assumed that k is given. Using standard exponential search, 

22 this requirement can be removed. More precisely, for / = 1 log log n, set 

23 ki = 2^ . Run the above algorithm for = fco, fci, . . . . If the algorithm succeeds, 

24 report the result. If not, abort as soon as it turns out that an intermediate onion 

25 has more than ki layers and try A'i+i. The total time is 

log log k 

0{n2') ^0{nlogk), 

1=0 

26 as desired. This finally proves our main result. 

27 Theorem 3.4. Let V be a set of n disjoint unit disks in M^. We can build a 

28 data structure that stores V, of size 0{n), in 0{n\ogn) expected time, such that 

29 given a sample P ofV, we can compute ^{P) in 0(n\ogk) time, where k is the 

30 number of layers in ^{P). □ 

31 Remark. Using the same approach, without the exponential search, we can 

32 also compute the outermost k layers of an onion with arbitrarily many layers in 

33 0{n\ogk) time, for any k. In order to achieve this, we simply abort the union 



34 algorithm whenever k layers have been found, and note that by Observation 3.1 

35 the points in P not on the outermost k layers of '${P) will never be part of the 

36 outermost k layers of %i{Q) for any Q D P. 
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1 4 Deterministic Preprocessing 



2 We now present alternatives to Lemma |3.4| First, we describe a very simple 

3 construction that gives a deterministic way to build an (9/10 + e, 1/2 + e)-SDT 

4 in 0{n\ogn) time. 

5 Lemma 4.1. Let V be a set of m non-overlapping unit disks. Suppose that the 

6 centers of T> have been sorted in horizontal and vertical direction. Then we can 

7 find in 0{m) time a (vertical or horizontal) line ('., such that f intersects 0{\/m) 

8 disks and such that (' has at least m/l{) disks from T) completely to each side. 

9 Proof. Let Vi, V,., Vt, V,. be the m/lO left-, right-, top-, and bottommost disks 

10 in V, respectively. We can find these disks in 0{rn} time, since we know the 

11 horizontal and vertical order of their centers. We call Vq :— 2?; U Vr WDtWDr 

12 the outer disks, and 2?,; := V \ the inner disks. 

13 Let R be the smallest axis-aligned rectangle that contains all inner disks. 

14 Again, R can be found in linear time. There are f](m) inner disks, and all disks 

15 are disjoint, so the area of R must be i7(m). Thus, R has width or height f2{y/Tn); 

16 assume wlog that it has width f2{^m). Let R' C Rhe the rectangle obtained by 

17 moving the left boundary of R to the right by two units, and the right boundary 

18 of R to the left by two units. The rectangle R' still has width f2{s/m), and it 

19 intersects no disks from Vi U V,.. There are Q{-\/m) vertical lines that intersect 

20 R' and that are spaced at least one unit apart. Each such line has at least m/10 

21 disks completely to each side, and each disk is intersected by at most one line. 

22 Hence, there must be a line that intersects 0{y/m) disks, as claimed. We can 

23 find such a line in 0{m) time by sweeping the disks from left to right. □ 

24 The next lemma improves the constants of the previous construction. It 

25 allows us to compute an (1/2 + e, 5/6 + t)-SDT tree in deterministic time 

26 O(nlog^n), but it requires comparatively heavy machinery. 

27 Lemma 4.2. Let V be a set of m congruent non- overlapping disks. In determin- 

28 istic time O (to log to), we can find a line ( such that there are at least m/2~cm^/^ 

29 disks completely to each side of £. 

30 Proof. Let X be a planar n-point set, and let 1 < r < n be a parameter. A 

31 simplicial r -partition of X is a sequence zii , . . . , Z\r of r triangles and a partition 

32 X = Xi(j ■ ■ ■ (jXr of X into r pieces such that (i) for « = 1, . . . , r, we have Xi C 

33 Ai and \Xi\ e {n/r, . . . ,2n/r}; and (ii) every line {' intersects 0{y/r} triangles 

34 Ai. Matousek showed that a simplicial r-partition exists for every planar n-point 

35 set and for every r. Furthermore, this par titio n can be found in O(nlogT) time 



36 (provided that r < n^^^ , for some (5 > 0) 18 Theorem 4.7]. 

37 Let a, /? > be two constants to be determined later. Set r := to". Let Q be 

38 the set of centers of the disks in V. We compute a simplicial r-partition for Q in 

39 0{m log m) time. Let Z\i, . . . , Z\„ be the resulting triangles and Q = QiU ■ ■ • OQr 

40 the partition of Q. Set s := , and for i = L . . . , .s, let be the line through 

41 the origin that forms an angle (i/2s)7r with the positive x-axis. Let Yi be the 
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1 projection of the triangles Ai, . . . , A,, onto . We interpret Y; as a set of weighted 

2 intervals, where the weight of an interval is the size \Qj\ of the associated point 

3 set for the corresponding triangle. By the properties of the simplicial partition, 

4 the interval set Y) has depth 0{^), i.e., every point on is covered by at most 

5 0{\/r) intervals of Y,. 

6 Note that the sets Y, can be determined in 0(.srlogr) = 0(m"+'^ log m) = 

7 0{m) total time, for q, d small enough. Now, for each Y;, we find a point ci on 

8 that has intervals of total weight m/2 - 0{y/r{m/r)) = m/2 - ©(m^-^/^) 

9 completely to each side. Since the depth of Yj is 0{y/r), we can find such a point 

10 in time O(logr) with binary search, for a total of O(.slogr) = 0{in) time (it 

11 would even be permissible to spend time 0(r) on each 1";). Let be the line 

12 perpendicular to f' through r,. 

13 The analysis of Alon et al. shows that for each , there are at most 0{s log .s) 

14 disks that intersect and at least one other line £j [2j Section 2]. Thus, it 

15 suffices to focus on the disks in V that intersect at most one line £i. By simple 

16 counting, there is a line f, that exclusively intersects at most m/s — m^~^ 

17 disks. It remains to find such a line in 0{m) time. For this, we compute the 

18 arrangement A of the strips with width 2 centered around each together with 

19 an efficient point location structure. For each cell in the arrangement, we store 

20 whether it is covered by 0, 1, or more strips. Using standard techniques, the 

21 construction takes 0{s^) — 0{m'^'^) time. We locate for each triangle Ai the 

22 cells of A that contain the vertices of A^. This needs C>(r log. s) = 0(m" logm) 

23 steps. Since every line intersects at most 0{y/r) = 0{m°'/^) triangles, we know 

24 that there are at most 0{srn"^^) = 0{m^^°/^) triangles that intersect a cell 

25 boundary of A. We call these triangles the bad triangles. 

26 For all other triangles Ai, we know that the associated point set Q, lies 

27 completely in one cell of A. Let Vi be the corresponding disks. By using the 

28 information stored with the cells, we can now determine for each disk D G 2?; 

29 in 0(1) time whether D intersects exactly one line Thus, we can determine 

30 in total time 0{;m) for each line f, the total number of disks that intersect only 

31 £i and whose center is not associated with a bad triangle. Let £ be the line for 

32 which this number is minimum. 

33 In total, it has taken us 0{m\ogm) steps to find £. Let us bound the number 

34 of disks that intersect £. First, we know that there are at most 0(m'^+"/-^ • 

35 m^^") = 0(m^+'^^"/'^ ) disks whose centers lie in bad triangles. Then, there are 

36 at most 0(m'' log m) disks that intersect £ and at least one other line. Finally, 

37 there are at most m^^'^ disks with a center in a good triangle that intersect only 

38 £. Thus, if we choose, say, [3 — 1/6 and a — 2/3, then ( crosses at most 0{m^^^) 

39 disks in V. Furthermore, by construction, C has at least m/2 — 0{m^^'^) disk 

40 centers on each side. The result follows. □ 



41 Remark. Actually, we can use the approach from Lemma |4.2| to compute an 

42 (1/2 + e, 5/6 + e)-SDT in total deterministic time O(nlogn). The bottleneck lies 

43 in finding the simplicial partition for Q. All other steps take 0{m) time. However, 



44 when applying Lemma 4.2 recursively, we do not need to compute a simplicial 



45 partition from scratch. Instead, as in Matousek's paper, we can recursively refine 
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Fig. 4. The lower bound construction consists of (;/3 unit disks centered on a horizontal 
line (5 in the figure), and two groups of n/3 points sufficiently far to the left and to the 
right of the disks. Distances not to scale. 



Fig. 5. ii/k copies of the construction on a regular n/fc-gon. 



18 Corollary 3.5] (while duplicating the 



1 the existing partitions in hnear time 

2 triangles for the disks that are intersected by (). Thus, after spending 0{n\ogn) 

3 time on the simphcial partition for the root, we need only linear time per node 

4 to find the dividing lines, for a total of 0{n\ogn), by Lemma 3.3 



5 5 Lower Bounds 



6 We now show that our algorithm is optimal in the decision tree model. We begin 

7 with a lower bound of i7(nlogn) for k — n{n). Let n be a multiple of 3, and 

8 consider the lines 

C : y = -1/2 - 6/n - x/n^; C : y ^ -1/2 - 6/n + x/n^. 

9 Let Vn consist of ti/3 disks centered on the .i--axis at :r-coordinates between — n/6 

10 and n/6; a group of «/3 disks centered on £^ at .x-coordinates between and 

11 Ti^ + n/3; and a symmetric group of ri/3 disks centered on f+ at :r-coordinates 

12 between —n^ — n/3 and —n^. Figure E] shows Pis. 

13 Lemma 5.1. Let tt be a permutation on n/3 elements. There is a sample P of 

14 V„ such that pi (the point for the ith disk from the left in the main group) lies 

15 on layer n[i) of%{P). 
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1 Proof. Take P as the n/3 centers of the disks in V on f^, the n/3 centers 

2 of the disks in V on £+, and for each disk Di E V on the :r-axis the point 

3 Pi — {i — n/6,7r(j) • 3/n — 1/2). By construction, the outermost layer of 

4 contains at least the leftmost point on £+, the rightmost point on f^, and the 

5 highest point (with (/-coordinate 1/2). However, it does not contain any more 

6 points: the line segments connecting these three points have slope at most Ijv? . 
1 The second highest point lies 3/??, lower, and at most n/3 further to the left or 

8 the right. The lemma follows by induction. □ 

9 There are ('ri/3)! — 2'='("'°s") permutations tt; so any corresponding decision 

10 tree has height Q(n\o^n). We can strengthen the lower bound to i7(nlogfc) by 

11 taking njk copies of and placing them on the sides of a regular (n/fc)-gon, 

12 see Figure [5] By Lemma |5.1[ we can choose independently for each side of the 

13 (n/fc)-gon one of (fc/3)! permutations. The onion depth will be A'/3, and the 

14 number of permutations is ((/s/3)!)"/''' = 2^^("'°s''). 

15 Theorem 5.1. Let e N and n > k. There is a set V of n disjoint unit disks in 

16 E^, such that any decision-has ed algorithm to compute i${P) for a sample P of 

17 V, based only on prior knowledge of T) , takes Q{n\ogk) time in the worst case. 

18 The lower bound still applies if the input points come from an appropriate 



19 probability distribution (e.g., [l| Claim 2.2]). Thus, Yao's minimax principle 19 



20 Chapter 2.2] yields a corresponding lower bound for any randomized algorithm. 



21 6 Conclusion and Further Work 

22 It would be interesting how much the parameter k can vary for a set of imprecise 

23 bounds and how to estimate A: efhciently. Further work includes considering more 

24 general regions, such as overlapping disks, disks of different sizes, or fat regions. 

25 It would also be interesting to consider the problem in 3D. Three-dimensional 

26 onions are not well understood. The best general algorithm is due to Chan and 

27 needs O(nlog^n) expected time [H], giving more room for improvement. 
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